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Abstract 

A new geometry to trap neutral particles with an ac electric field using a simple electrodes structure is 
described. In this geometry, all electrodes are placed on a single chip plane, while particles are levitated 
above the chip. This provides an easy construction of the trap and a good optical access to the trap. 

particle storage, atom optics, laser cooling 

1 Introduction 

To investigate the properties of particles such as atoms and molecules, it is desirable to hold them in an 
isolated environment. For this purpose, various kinds of traps have been developed. Among them, magneto- 
optical trap (MOT) and magnetic trap (MT) are successfully used for a variety of applications [T]. However 
they require the particles to have either a closed transition (MOT) or a magnetic moment (MT), which 
restrict the range of particles to be trapped. Electric trap (i.e. a trap that uses electric field), on the other 
hand, has the advantage that it can trap almost any kind of neutral particles. It is easily shown that it is 
impossible to trap neutral particles using a static electric field [2], and an ac electric trap was proposed to 
overcome this difficulty. [31 21 [5] Disadvantage of such electric trap is its shallow trap depth (potential depth 
for sodium atoms is 190/iK for an electric field of lkV mm -1 ), and one needs to make a small trap to increase 
the curvature of the potential and hold particles against the gravity. However ac electric traps proposed so 
far had solid structures and the fabrication of a small trap was rather difficult. For such technical reasons, 
ac electric traps were not realized until recently. [6j [Tl [8j [9] 

In this paper, we propose a new geometry of an ac electric trap in which all electrodes are placed on 
a single plane, while particles are trapped above that plane. Such geometry enables one to fabricate the 
trap structure in mass with high precision, which opens the possibility of making micro-structured atom 
and molecule chips that benefit from the advantages of long decoherence time [10], low electric power 
consumption, and scalability of the electric traps. 

2 Principle of the dynamical trapping 

We start from the trap configuration described in [[5]], schematically shown in fig. [1] in which all electrodes 
sit on a single plane (xy-plane) . Two pairs of electrodes are placed on x- and j/-axis symmetrically (i.e. at 
(±so,0, 0) and (0, ±so,0)). There are two phases in the operation of this trap: In phase A (phase B) a 
voltage is applied on the pair of electrodes on a;-axis (j/-axis). The potential that a neutral particle feels 
under the electric field E(r) is 

l/(r) = -ia|E(r)| 2 (1) 

where a is the polarizability of the particle, a is always positive for atoms or molecules in a stable state. In 
phase A or B, the potential V(r) has the form 

V A (r) = \miul{~i l0 x 1 + y 2 + i z 2 ) + V 
V B {v) = §mwg(+x 2 -W , + £o* 2 ) + Vb 

up to 2nd order in r. Here, m is the mass of the particle and uiq is the angular frequency of oscillation in 
y-direction in phase A (a;-direction in phase B). For spherical electrodes r/o = 2 and £o = By switching 

between phase A and B alternatively for time Ta and Tb (usually Ta = Tb), particles are trapped around 
the origin (0, 0, 0): They are dynamically captured in xy-plane (similar to the RF ion trap [H]), whereas the 
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Figure 1: Two pairs of electrodes are alternatively switched on and off between phase A and B. [phase A] 
Force is attractive in yz-plane and repulsive in ^-direction, [phase B] Attractive in zx-plane and repulsive in 
y-direction. 
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Figure 2: Trap point is the saddle point of the density of electric field lines. By adding two additional charges, 
another saddle point is created above xy-plane (right). 



confinement in ^-direction is static. Larger £o gives stronger confinement in ^-direction, while the stable 
region for the driving frequency Q — 2n/T is reduced with increasing r/o (T = Ta + Tb is the period of the 
applied voltage). 



3 How to lift up the trap point from the chip plane 

Particles are trapped at the saddle point of j_E(r)j 2 . In order to lift up this trapping point off the :ry-plane, 
we place additional electrodes in between the existing electrodes to bend the field lines and make a saddle 
point other than the origin (fig. [2]). We still keep the electrodes configuration symmetric under x <-> —x, 
y <-> — y, x <-> y. Now (f2| becomes 



Va(t) = lmu a {-Tp?+v'+t(z-h) a } + V, 
V B (v) = \uluj 2 {+x 2 -w 2 +t{z-K) 2 } + V. 



(3) 



near the new saddle point (0, 0, h). Let <f>(r) be the scalar potential. In phase A (and similarily in phase B), 
because of 4>{x, y, z) = —cj>(—x,y,z), 4>{x, —y, z) = <f)(x,—y,z), 



E x (x,y,z) = +E x (-x,y,z) 
E y (x,y,z) = -Ey(-x,y,z) 
E z {x,y,z) = -E z (-x,y,z) 
E x (x,y,z) = +E x (x,-y,z) 
E y (x, y, z) = -E y (x, -y, z) 
E z (x,y,z) = +E z (x,-y,z) 



(4) 
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— and V x E = 0. From these equalities, 
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) =0, and thus rj — £ = 1 is still satisfied. 



3.1 Design procedure: step 1 

First, we consider electrodes as point charges and place them on xy-plane as shown in fig. [3] Four 
outer point charges placed at s ou t = (±so,0), (0, ±so) correspond to the four electrodes of the origi- 
nal configuration (fig[T]) and we put charges of ±go or at appropriate phase. Four inner point charges 
at Si n = (±s,±s), (±s,=ps) are added to the original configuration with charges ±g depending on the 
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Figure 3: In addition to four outer point charges that corresponds to the four electrodes of of the original trap 
configuration, we place four inner point charges. 
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Table 1: parameter comparison table for s = 0.3sq- 



phase (see figure). In fig. [4] we plot parameters £ that appear in ([3| as functions of ^- for sev- 

eral j- values (77 is calculated from ( as f/ = ( + 1). Here u>o is u for the original configuration, i.e. 

= y/^V A \ ^ = ^y^i^T- Normalized potential V norm (r) = 2 A (l^)V(r) is plotted 
along z-axis for different q in fig. [S] for different s in fig. [1J] and along x- (y-) direction in fig. [7] We choose 
j- = 0.3 and — 0.06 to proceed the design procedure further. Relevant parameters for j- — 0.3 and 

= 0.06 are compared with those for the original configuration in tableQ] (see section^for x). 

3.2 Design procedure: step 2 

Now we are going to replace the point charges by electrodes of finite size. For that, first we calculate 
the equipotential surface of the electric field. Figure [8^, is the plot of the normalized scalar potential 
4>norm(r) = ^gp" <ft( r ) m the chip plane (z — 0). Equipotential surfaces of (finorm = ±a (a > 0) around the 
outer point charges at s ou t are nearly spherical for large a and small q/qo, and can be well imitated by a 
disc shaped electrodes of radius r ou t = ^a+l cen tered at s ou t- However smaller a gives higher electric field 
assuming that the voltage applied on the electrodes is fixed. In fig. [5b we choose a = 2.0 (r out = 0.4so). 
We apply to these outer electrodes voltage of ±«o or at appropriate phase: 4jTC ° s ° = f- (i.e. V{r) = 

20^ (so) 2 ^" orm ( r ))- Equipotential surfaces of <\> norra = around the inner point charges at Si n are again 
nearly spherical for ^ <C , and the inner point charges can be replaced by disc shaped electrodes of radius 

r in = ^-{(sq - 2sos + 2s 2 )~i - (sq + 2s s + 2s 2 )~i}~ 1 (r ln = 0.11s for q = 0.06g , s = 0.3s ) centered at 
Sin which are always connected to the ground (0 voltage). 



4 Trapping particles under the gravity of the earth 

In this section, we give some practical values to trap neutral atoms and dielectric spheres on the earth. We 
orient z-axis along the gravity direction. We set a condition that the gravitational sag should be smaller 
than as (with a < 1): < as . This gives uj > J as ^j^. ■ On the other hand u = ji\f^^- We 

fix here the trap design as = 0.3, ^ = 0.06 and a = 2.0 which gives ^ = 0.85, £ = 0.94. We also put 
a = 0.1 and g = 9.81 m s~ 2 . 

First we discuss on trapping neutral atoms. We choose sodium atoms as an example: a = 2.68T0 _39 F m 2 , 
[T5] m = 3.82 x 10" 26 kg, and thus ^ = 7.0 x 10~ 14 Fm 2 kg~\ If we set s = 0.5mm, then the condition on 
the gravitational sag requires ujo > 2n x 79Hz. This will be satisfied with the applied voltage of Vo > 270V. 
Setting vo = 270V, then u> = 2n x 67Hz and, using fig. 1121 the trap should be stable for 7.8ms < T < 9.5ms 
with T A =T B = f . 

Next we discuss on trapping of a dielectric material. The polarizability of a dielectric sphere of radius a 
and dielectric constant e is a = 47ren "T^ a 3 (mass is m — ^-p a 3 where p is the density) so that — = ^£fl. e ~ e Q 
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Figure 4: Plot of the trap parameters: relative trap height ^, relative trapping strength ij, and £ are plotted 
as functions of -2- for =0.2, 0.3, and 0.4. 
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Figure 5: Plot of the Stark potential along z-axis for different q. 




Figure 6: Plot of the Stark potential along z-axis for different s. 




Figure 7: Plot of the Stark potential along x-direction for phase A and phase B (s = 0.3s , q = 0.06qo). 
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Figure 8: (a) Cross section of the cquipotcntial surface in phase A: normalized scalar potential 0„ O rm(r) for 
s = 0.3s an d q = 0.06qo at the chip surface (z = 0) is plotted as a function of x and y. (b) Disc shaped 
electrodes shown in black imitate the equipotential surfaces. Inner electrodes are connected to the ground. 
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Figure 9: Plot of y as a function of — for several — values. 

is independent of its size. Polystyrene(e = 2.5eo, p = 1-05 x 10 3 kgm -3 ), for example, has the value 
^ = 8.4 x 10 _15 Fm 2 kg _1 , which is roughly one order smaller than that of sodium atoms. 

5 Trapping of polar particles 

In the case of polar particles of permanent dipole moment pL, the force under the electric field can be written 
as, assuming the direction of the dipole moment is always oriented parallel to the electric field (see appendix 

fi(r)= v. /',<•;,/•:.;.•, 

- mE^^W (5) 

Thus the force is derived from the potential 

V(r) = - M £(r) (6) 

Near the trap center r c , E(v) = E< ~ r ^ 2 E(r^)' : ^ > so ^ aa ^ the dynamics of the polar particles can be argued in the 
same way as that of the nonpolar particles by replacing a with E ^ r ? . We define a dimensionless parameter 

X ee where E ee E(Q)\ q=0 = ^^IV^ (0)1^=0 = ^- Now u = ^^St* = 

A plot of \ is given in fig[9] For a particle of bulk material of permanent polarization P and density p, 

^ = ^ so that Wo = \J~^~f§z i s independent of its shape and size. To give practical parameters, we again 

use the trap design s = 0.3so, q = 0.06go, and a = 2.0, and consider trapping of BaTiOs micro particles 
(P = 0.26Cm" 2 and p = 5.5 x 10 3 kgm~ 3 ). If we set so = 3.0mm, then the condition on the gravitational 
sag with a = 0.1 requires ujq > 2n x 32Hz which implies vo > 5.0V. At vo = 5.0V, lo = 2-n x 27Hz and, from 
fig. 1121 the trap should be stable for 19ms < T < 23ms. 

6 Conclusions and outlook 

We have presented a new design of ac electric trap for neutral particles and have shown that it is possible 
to make a trapping point above the chip surface with planar electrodes structure. We have also discussed 
about the feasibility of this design by showing typical parameters for trapping neutral atoms, dielectric 
spheres, and polar particles. The parameters for this trap were comparable with those demonstrated in [[6]]. 
In the final step of the design, we have used disc shaped electrodes to imitate the equipotential surface of 
the electric field. [14] To analyse its consequence, we are now calculating the electric field using numerical 
methods. The result is still preliminary, but it has turned out that the field near the trapping point is not 
so much dependent on the shape of the inner electrodes. Thus, as an application of our trap, a particle 
conveyor such as shown in fig. [TO] should be also possible. 
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Figure 10: Possible design of a particle conveyor. '+' marks show the trapping points levitated off the chip 
surface. Electrodes shown in black are connected to the ground. 



phase 




Figure 11: Time evolution matrix U. 
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A Stability 

In this section, we investigate the stable condition of the dynamical confinement. The equations of motion are 
separated in x-, y-, ^-directions so that we explore the motion in ^-direction only. We follow the procedure 
described in [[4]] and dehne a state vector X(t) as 

= ( $j ) « 

where x(t) (v(t)) is the position (velocity) of the particle at time t. Time evolution of X is described by a 
time evolution matrix U(t): X(t) = U(t)X(0). 

a?W(t) x {2 \t) 

'(*) « (2) w 



with 17(0) = I (I is the unit matrix). We assume that the phase A starts at t = (hg. 1 1 1 p so that 

U(nT + t) = U(r)U(T) n with U(T) = UbUa where Ua (U B ) is the time evolution matrix for the phase A 
(B)- 

Ua = ( cosh^/TjujTA ^sinh^TA \ 

\ y/rju sinh y/rjujTA cosh y/rjujTA J 

( cos^Tb ^ sin uTb \ (1Q) 

V — cjsmuj 1b cos luIb J 

with Ta — Tb — . Let A be the eigenvalue of U(T). 



X 2 — 2,3 A + 1 = 
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Figure 12: Plot of j3 as a function of ojT and rj. The trap is stable in the region where \j3\ < 1. 



with /3 = ±Tr U(T) = cosh ^r\ujT A cos wTs + sinh ^luT a sin wT B . From (TT} 

iAi-1 1 (l/3i - 1} ri2i 

|A| -\ \0\±^i m>i) (12) 

The trap is stable if and only if |/3| < 1. In fig. [12] we plot /3 as a function of tuT and 77 (T A = T B = f ) to 
show the stable region. 



A.l Unharmonicity 

So far we have not included the unharmonic terms in the potential. In the original Stark chip configuration, 
the lowest unharmonic term was 4th order because of the z — > — z symmetry. In our configuration, due to 
the absence of this symmetry, 3rd order terms are present. However these terms have the form 

V 3 (x,y,z) = (71 x 2 + 7 2 t/ 2 )(z - h) +yi{z~hf (13) 

and are still harmonic in the plane z = const, (i.e. the plane in which the confinement is dynamical), we 
expect that the 3rd order terms do not reduce the stable region in phase space drastically. 



B Trap tightness 

To estimate the tightness of the dynamical confinement, we apply a constant force Fo in ^-direction on 
the particle in the trap and calculate the shift of the trap center (i.e. the time average of the position of 
the particle). Let W{i) be the state vector at time t with the initial condition W(0) = 0. Then the time 
evolution of a general state vector is now X(t) — U(t)X(0) + W(t). 

W(2T) = U(T)W(T) + W(T) 

W(3T) =U(T){U(T)W(T) + W(T)} + W(T) 

W(nT) ={U{T) n - 1 + U(T) n - 2 + ... + U(T) + l}W(T) ( 14 ) 
= ¥^§f W(T) 

1 - U(T) n 

W(nT + T ) = U(t) ] _ _ J ' W(T) + W(r) (15) 

X(nT + r) = U(T)U(T) n X(Q) + U(t)^^W(T) + W(r) (16) 

To evaluate the time average of X(i), we put t = nT + r with n = 0, 1, 2, ... and < r < T, and take the 
average over n and r. Because U(T) n = inside the stable region, we obtain 



x(nT + r) 



X(nT + r) = ^ -^-L j = U(r) Y - m W { T) + W(r) (17) 

After a lengthy calculation, fllT|) yields 

x(nT + r)^ { -y + v^^-i) s in^-^7(co^r B -i) sin h ^t a (i + !) 2 | (18) 
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Figure 13: Plot of ^ as a function of wT and r\ (T A =T B = f ). 
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C Dynamics of the polarization orientation of a polar parti- 
cle 

Consider a polar particle of size L and volume u (u ~ L 3 ) that consists from a bulk material under the 
electric field E. Let 6 be the angle between the electric field E and the polarization P of the material. Then 



N 



-uPE sin 



(19) 



where / is the moment of inertia and iV is the torque. I ~ pL where p is density of the material (/ = f^pL 5 



for a sphere of diameter L) . Assuming 9 is small, 9 
oscillation of the polarization direction is thus u) osc 
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"pi 1 
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L vu± ,^v>. ~ osc — , , — r ). If the switching 
time between phase A and B is longer than aj~/ 3 c> then the direction of the polarization adiabatically follows 



that of the electric field. For BaTiC>3 (see section [5| with L 
(sphere: lo osc = 7.8 x 10 5 s~ 1 ). 
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